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ABSTRACT 

From the Sloan Digital Sky Survey (SDSS) Data Release 12, which covers the full 
Baryonic Oscillation Spectroscopic Survey (BOSS) footprint, we investigate the possi¬ 
ble variation of the fine-structure constant over cosmological time-scales. We analyse the 
largest quasar sample considered so far in the literature, which contains 13 175 spectra 
(10 363 from SDSS-III/BOSS DR12 + 2812 from SDSS-II DR7) with redshift z < 1. We 
apply the emission-line method on the [O hi] doublet (AA 4960,5008 A) and obtain A a/a = 
(0.9 ±1.8) x 10 3 for the relative variation of the fine-structure constant. We also investi¬ 
gate the possible sources of systematics: misidentification of the lines, sky OH lines, H fi and 
broad line contamination, Gaussian and Voigt fitting profiles, optimal wavelength range for the 
Gaussian fits, chosen polynomial order for the continuum spectrum, signal-to-noise ratio and 
good quality of the fits. The uncertainty of the measurement is dominated by the sky subtrac¬ 
tion. The results presented in this work, being systematics limited, have sufficient statistics to 
constrain robustly the variation of the fine-structure constant in redshift bins (Az ~ 0.06) over 
the last 7.9 Gyr. In addition, we study the [Ne hi] doublet (AA 3869,3968 A) present in 462 
quasar spectra and discuss the systematic effects on using these emission lines to constrain the 
fine-structure constant variation. Better constraints on A a/a (< 10“ 6 ) using the emission-line 
method would be possible with high-resolution spectroscopy and large galaxy/qso surveys. 

Key words: line: profiles - quasars: emission lines - cosmology: observations - surveys - 
large-scale structure of Universe. 


1 INTRODUCTION 

Since Dirac’s philosophical argument (Dirac 1937) against the 
fixed value of fundamental constants of Nature, several experi¬ 
ments have been performed to constrain possible variation on di- 
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mensionless constants of physical theories. Fundamental constants 
of physics could be thought of as parameters which enter in our 
description of Nature but they cannot be predicted with our cur¬ 
rent theories and should be measured. Dirac’s idea is based on the 
unlikely fact that the most fundamental constants of the Universe 
have a certain fixed value (at a given energy) with no apparent re¬ 
lation with the real world. It is more likely that their present val¬ 
ues are the result of a dynamical process, which had yielded the 
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Table 1. Summary of the results obtained by recent works based on the [O hi] emission line method for the possible variations of the line-structure constant. 


Reference 

Quasar spectra 

SDSS release 

£min 

Zmax 

Time ago (Gyr)* 11 ' 

A a/a (xl0~ 5 ) 

Bahcall et al. (2004) 

42 

EDR (Stoughton et al. 2002) 

0.16 

0.80 

7.0 

7 ± 14 

Gutierrez & Lopez-Corredoira (2010) 

1568 

DR6 (Adelman-McCarthy et al. 2008) 

0.00 

0.80 

7.0 

2.4 ±2.5 

Rahmani et al. (2014) 

2347 

DR7 (Abazajian et al. 2009) 

0.02 

0.74 

6.7 

-2.1 ± 1.6 

This work (2015) 

13 175 

DR12 (Alam et al. 2015) 

0.04 

1.00 

7.9 

0.9 ±1.8 


For a ACDM cosmology with Hq = 67.8 kms -1 Mpc -1 , Cl m = 0.31 and Q.\ = 0.69 from Planck+WMAP- 9+BAO (Planck Collaboration et al. 2014). 

® Note: Since we have a larger sample than Gutierrez & Lopez-Corredoira (2010), we expect a factor k 2.5 of improvement in the error just from purely 
statistical reasons. In Figs 9 and 10, it is shown that the error is dominated by the sky subtraction algorithm, which suggests that the performed analysis have 
reached the maximum precision with the available data. 


fundamental constants as they are measured today. Therefore, they 
should be considered as characterizing the state of the Universe 
(Uzan 2003). There are many current theoretical frameworks which 
allow for such variation of the fundamental constants, for instance, 
string theory (Maeda 1988), modified gravity and theories with ex¬ 
tradimensions (e.g. Clifton et al. 2012). Moreover, the experimental 
bounds on their variation have become a stringent test for those the¬ 
oretical models (e.g. Thompson 2012; Leal et al. 2014). The most 
studied fundamental constants are the fine-structure constant a , the 
Newton gravitational constant G and the electron-to-proton mass 
ratio (Uzan 2003, 2011; Garcla-Berro et al. 2007). 

The fine-structure constant governs the electromagnetic cou¬ 
pling between photons and charged particles a = e 2 /(hc). The cur¬ 
rent constraint on its relative variation A a/a, over geological time- 
scales, is |Ao-/a| < 7 x 10 -8 up to z » 0.15 (2Gyr ago); obtained 
from the Oklo phenomenon (e.g. Petrov & et al. 2006). It has also 
been reported |Aff/a| < 3 x 10 -7 up to z ~ 0.45 (4-5 Gyr ago) from 
meteorites (Olive et al. 2002); which also excludes possible varia¬ 
tions on the scales of the Solar system. On the other hand, there are 
also constraints, |Aor/or| < 10 2 , based on the cosmic microwave 
background (CMB; Landau & Scoccola 2010; Planck Collabora¬ 
tion et al. 2014) at z ~ 1100 and from big bang nucleosynthesis, 
the latter being model-dependent. By measuring fine-structure mul- 
tiplets at different redshift in the absorption or emission spectra of 
galaxies and quasars, located at different directions in the sky, one 
can measure an estimate of the variation of a with time or space 
over cosmological scales. 

The first measurements on the variation of a from astronom¬ 
ical observations reached an accuracy of A a/a ss 1CT 2 - 1(L 3 
(Savedoff 1956; Bahcall & Salpeter 1965; Bahcall & Schmidt 1967; 
Bahcall, Sargent & Schmidt 1967). Since then, the methodology 
and understanding of systematics has dramatically improved. Cur¬ 
rent measurements of absorption multiplets along the line of sight 
of three quasars around redshift 1.5, observed with spectral resolv¬ 
ing power R » 60 000 at UVES/ESO-VLT, reached the ss 5 x 10 -6 
level (Evans et al. 2014). Using emission lines, an accuracy of 
2 x 10 -5 was achieved analysing 1500 - 2300 quasar spectra 
at z ~ 0.6 (Gutierrez & Lopez-Corredoira 2010; Rahmani et al. 
2014), taken with the Sloan Digital Sky Survey (SDSS) R ~ 2000 
spectrograph. 

The measurements on absorption features on a quasar spec¬ 
trum are currently limited by the precision in the absolute wave¬ 
length calibration of the spectra, i.e., 50 - 200 ms -1 using spectra 
with R 60,000 (Molaro et al. 2013; Evans et al. 2014; Whit¬ 
more & Murphy 2015). Furthermore, the so-called many-multiplet 
(MM) method used in Evans et al. (2014), although more precise, 
remains controversial as several assumptions are made, the most 


important one being ionization and chemical homogeneity. These 
assumptions may induce systematic biases on the value of a. 

In this article, we use the method based on the [O in] emission 
lines, first proposed by Bahcall & Salpeter (1965), which is less af¬ 
fected by systematics. In particular, there is no need for assuming 
ionization and chemical homogeneity, since the studied lines have 
the same profile (the transitions originate at the same upper energy 
level). Furthermore, the emission-line method suffers of much less 
spectral distortion, since the measurements of A a/a are done on 
a spectral window ~ 100 A as compared to ~ 1000 A when the 
MM method is used. With a large ensemble of quasars and/or us¬ 
ing high-resolution spectroscopy, the uncertainty can be reduced 
significantly, and will compete with the absorption method when 
using high-resolution spectroscopy. 

The beginning of the SDSS survey opened a new era of pre¬ 
cision, allowing us to use big samples of quasars; thus, reducing 
the statistical uncertainty of the measurement of A a/a (see Ta¬ 
ble 1). Here, we extend these works by using the SDSS-III/BOSS 
Data Relase 12 (SDSS-DR12; Alam et al. 2015), which covers 
the full Baryonic Oscillation Spectroscopic Survey (BOSS) sur¬ 
vey footprint with an area coverage of 10000 deg 2 . In contrast to 
these previous investigations, we use spectra obtained with the cur¬ 
rent BOSS spectrograph (Smee et al. 2013) instead of the previ¬ 
ous SDSS-I/II instrument, making our BOSS sample totally inde¬ 
pendent from previous works. Moreover, the spectral range of the 
BOSS spectrograph allows an extension of the redshift interval for 
the [O in] doublet from z = 0.8 to z = 1. The number of quasar 
spectra is increased by a factor of 5 with respect to SDSS-DR7. All 
these spectra have been visually inspected and classified as quasars 
by the BOSS collaboration, and their products are provided in the 
SDSS-III/BOSS Data Release 12 Quasar catalogue (DR12Q; see 
Paris et al. 2015). For the final constraint on A a/a, we combine in 
this work the BOSS sample with the previously studied SDSS-II 
DR7 quasar sample. 

There are several emission doublets, in addition to [O hi] 
(AA 4960,5008 A), that can be used to measure A a/a as noted by 
Bahcall et al. (2004), and first used by Grupe et al. (2005). Gutierrez 
& Lopez-Corredoira (2010) analysed different doublets and found 
that the [Ne hi] (AA 3869,3968 A) and [Sin] (AA 6719,6733 A) 
doublets appear in quasar spectra with sufficient frequency to have 
a meaningful sample. Results for [Si n] are consistent with no vari¬ 
ation of the fine-structure constant, although the uncertainty is an 
order of magnitude bigger than for [O hi] , and this doublet can 
only be used at low redshift < 0.4 for optical spectra. However, 
they obtained a positive variation of the fine-structure constant, 
A a/a = (34 ± 1) x 10~ 4 , when the [Ne m] lines are used. No expla¬ 
nation was found for this positive variation. In this work, we also 
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Figure 1. Left-hand panel: sky distribution of the full SDSS-III/BOSS DR12Q quasars (297 301) in J2000 equatorial coordinates. Right-hand panel: number 
of quasars with [O hi] emission lines in our fiducial sample (10363 quasars) in A z = 0.05 bins. S/NrQ m | 5008 > 10 (10363 quasars), black solid line; 
S/N[o in] 5008 > 25 (4015 quasars), blue dashed line; and S/Nj-q nl i 5008 > 50 (1498 quasars), red dotted line. 


analyse the [Ne hi] lines to check whether the same effect is present 
in our BOSS quasar sample. 

There are investigations which use Si rv absorption lines 
(/hi 1394, 1403 A) to obtain a precision of 4 x 10~ 6 (Chand et al. 
2005). This method also avoids the assumption of ionization and 
chemical homogeneity. However, since the separation between both 
lines is only *9A, the wavelength precision needed in the labora¬ 
tory for the separation between both lines is five times higher than 
using [O m] lines. Nevertheless, these constraints apply to the red- 
shift interval 1.59 < z < 2.92, which does not overlap with our 
range, thus they are complementary to the ones reported in this re¬ 
search. 

Finally, in the light of the upcoming large galaxy surveys, like 
eBOSS and DESI, that will provide millions of high-redshift galaxy 
spectra, we also discuss using galaxies instead of quasars to set 
constraints on the fine-structure constant. 

The paper is organized as follows. First, in Section 2, we de¬ 
scribe the data set used for our analysis. Next, in Section 3, the 
methodology is presented, the emission-line method is explained, 
and the code and simulations to analyse the spectra are described. 
In Section 4, we study several samples to check for systematics. 
Then, our results are presented in Section 5. Finally, we provide in 
Section 6 a summary of the main conclusions achieved with this 
research project. 


2 SAMPLE DESCRIPTION 

All the spectra used in this investigation were downloaded from 
the SDSS Database. This survey (York et al. 2000), which began 
taking observations in 1998, consists of a massive collection of op¬ 
tical images and spectra from astronomical objects including stars, 
galaxies and quasars. For this purpose, there is a dedicated 2.5- 
m wide-angle optical telescope at Apache Point Observatory in 
New Mexico (USA; for more details, see Gunn et al. 2006). The 
third phase of this project (SDSS-III; Eisenstein et al. 2011) in¬ 
cludes BOSS (Dawson et al. 2013) among its four main surveys. 
The data analysed in this research were provided by BOSS, and it 
is used for measuring A a/a for the first time. The SDSS-III/BOSS 
pipeline (Bolton et al. 2012) classifies the objects as quasars with 
&X 2 minimization procedure to fit the observed spectrum to multi¬ 


ple galaxy and quasar spectrum templates for all allowed redshifts. 
Then, a visually-inspected quasar catalogue is built from these ob¬ 
jects. Our fiducial sample is obtained from the DR12Q catalogue 
version (Paris et al. 2015). 

The wavelength coverage of the SDSS-III/BOSS spectrograph 
is 3600-10400 A and that of the SDSS-II spectrograph is 3800- 
9200 A. The BOSS sample is homogeneous since all the spectra 
have been obtained with the same instrument, and it is indepen¬ 
dent from previous investigations. The wider coverage of the new 
spectra allows consideration of higher redshifts (up to z = 1 for 
[O hi] doublet) than in the previous SDSS-II analysis based on the 
same method (see Table 1). The BOSS spectrograph has two chan¬ 
nels (blue and red) whose wavelength coverage is 3600-6350 A 
and 5650-10400 A, respectively. The resolving power ranges from 
1560 at 3700 A to 2270 at 6000 A (blue channel) and from 1850 at 
6000 A to 2650 at 9000 A (red channel). For our sample, the [O hi] 
lines fall in the red channel for > 96% of the quasars. The num¬ 
ber of pixels of each spectrum is about 4600 for the BOSS spectra 
and 3800 for the SDSS-I/II spectra. The pixel spacing is uniform in 
log-wavelengths (AlogT = 10~ 4 dex). More complete information 
about the SDSS-I/II and BOSS spectrographs can be found in Smee 
et al. (2013). 


2.1 Data selection 

The SDSS-III/BOSS DR12Q catalogue contains 297 301 objects. 
Fig. 1 (left-hand panel) shows the quasar distribution in the sky. 
We summarize below the main selection criteria in order to define 
our fiducial sample from this catalogue. 

(i) Redshift < 1. This limitation is imposed by the wavelength 
range of the BOSS optical spectrograph and the position of the 
[O hi] lines. This criterion decreases the sample down to 45 802 
quasars. 

(ii) S/Nj-q m ] 5008 > 10. We impose a mild constraint on the 
signal-to-noise ratio (S/N) of the stronger [O m] line (5008 A) in 
order to preserve a large number of spectra. Constraints on the ex¬ 
pected width and amplitudes of the lines help in avoiding misiden- 
tifications of the [O iii] doublet (see Section 4). This selection re¬ 
duces the sample from 45 802 to 13 023 objects. 






















4 


F. D. Albareti el al. 


A (A)(observed) 

6000 7000 8000 9000 



Figure 2. Left-hand panel: data points (1416) for which |e| = |5A Z /(1 + z) — <5Aol, namely the absolute value of the difference between the measured line 
separation at redshift z in rest frame and the local one, is bigger than 1 A plotted as a function of redshift (and the wavelength observed for [O iii] 4960). We 
compare with a typical sky spectrum: the [O iii] positions for these spectra correlate with the sky emission lines. Hence, these high values of |e| are due to bad 
sky subtractions and/or low S/N. These spectra are removed. Right-hand panel: a spectrum removed from the sample because of the sky emission-line criteria. 
For this quasar, we get e = 1.2 ± 0.6 A. The weak [O iii] line is affected by the subtraction of the 7995 A OH sky emission line, indicated by the vertical red 
dashed line. 


(iii) Non-converging fits. Since we analyse spectra with low 
S/N, there are some cases where the Gaussian fit to the lines does 
not converge. 1 244 spectra are discarded, leaving us with 11 779 
spectra. 

(iv) Sky emission lines. Strong atmospheric lines, for instance 
the O i 5578 A line, are poorly or not completely removed by the 
SDSS sky subtraction algorithm. This may lead to a wrong iden¬ 
tification of the [O iii] lines and to include low S/Nj-q m j spectra 
(Gutierrez & Lopez-Corredoira 2010). Both effects will produce 
outliers. We use the SDSS sky mask for Lycr forest studies which 
contains 872 lines (see Delubac et al. 2015, for more details) to 
remove spectra whose [O iii] lines lie within a particular distance 
from the strongest sky lines. Even though we vary the distance 
[O iii] - sky lines, use different set of sky lines (according to their 
intensity), or evaluate other conditions (S/N, fit errors, etc.) to re¬ 
move affected spectra; we usually eliminate 3-5 good spectra for 
each bad spectra eliminated. Thus, these tests decrease significantly 
the number of quasars while not being very effective: typically 50% 
of the outliers are not removed. Thus, we decided to eliminate all 
spectra for which the separation between both lines differ by more 
than 1 A from the local value (see the last paragraph in Section 3.3). 
Fig. 2 (left-hand panel) shows that the distribution of these outliers 
is correlated with a typical sky spectrum. From a visual inspec¬ 
tion, we observed that these spectra have low S/N, and they are 
in fact contaminated by sky emission line subtraction (see right- 
hand panel of Fig. 2). This effect causes us to discard 1 416 spec¬ 
tra (12% of the previous 11 779 quasars). Finally, we have 10363 
quasar spectra (our ‘fiducial sample’). 

The presence of broad H/3 emission line (4861 A) near the 
weak [O iii] line 4960 A could produce a blueshift in the deter¬ 
mination of the [O m] line position. This could mimic a positive 
variation on the fine-structure constant. Therefore, a constraint on 
the strength and/or width of the H /3 emission line has been im¬ 
posed on previous investigations (Bahcall et al. 2004; Gutierrez & 
Lopez-Corredoira 2010; Rahmani et al. 2014). However, we do not 
restrict any characteristic of the H/3 line in our fiducial sample. 


We obtain a weighted mean for Aa/a using as weights the uncer¬ 
tainty in Aa/a computed with the standard errors for the position of 
the lines derived from the Gaussian fits. The contamination of H/3 
is automatically taken into account. For instance, a broad H/3 line 
near the [O iii] 4960 line means a bad Gaussian fit. Thus, we obtain 
larger errors in the position of the line centroids and, consequently, 
in Aa/a. In Section 4, we analyse several samples where the S/Nh^ 
is constrained to check that the H/3 contamination has little weight 
on the final constraint value. 

An electronic table is published along with the paper which 
contains all the information of each spectrum from our fiducial 
sample of 10 363 quasars (see Appendix A). 

The distribution of the selected quasars in redshift accord¬ 
ing to their selected S/Njq m j 500g is plotted in Fig. 1 (right-hand 
panel). Fig. 3 (left-hand panel) displays a composite image built 
with all the spectra from our fiducial sample sorted by redshift. The 
right-hand panel shows the [O iii] doublet in rest frame. 

3 METHODOLOGY 
3.1 Measurement method 

To first order, the difference between the energy levels of an atom 
is proportional to a 1 . Transitions between energy levels of the same 
atom at a given ionization level, with the same principal quantum 
number and different total angular momentum J, have an energy 
difference proportional to a 4 . These groups of transitions are called 
fine-structure multiplets. Savedoff (1956) first realized that the fine 
structure of these energy levels could be used to break the degener¬ 
acy between the redshift effect and a possible variation of a. 

The value of the fine-structure constant can be measured 
through the separation between absorption or emission multiplets 
in the spectra of distant quasars (Uzan 2003) as 

Arr _ 1 | [(T 2 - Ti) / (T? + /li)] z f 
a 2 \ [(T 2 - Ti) / (T 2 + Ti)] 0 J 


( 1 ) 
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Figure 3. Composite image with our fiducial sample of 10 363 BOSS quasar spectra sorted by redshift. Left-hand panel: the whole range of wavelengths is 
shown. From right to left, the strongest emission lines are H a 6565 A; [O hi] AA 4960,5008 A; H/? 4861 A; H y 4341 A; [Ne iii] AA 3869,3968 A; [O n] 3730 
A; [Ne v] 3426 A; Mg n 2796 A and C iii] 1906 A. The narrow straight line at 5579 A is the strong [O i] atmospheric line. Right-hand panel: wavelength 
interval centred at the [O iii] doublet in rest frame. 


where A l2 (A 2 > Ai) are the wavelengths of the transitions and 
subscript 0 and z stand for their value at redshift zero (theoreti¬ 
cal/laboratory values) and at redshift z, respectively. For illustrative 
purposes, expression (1) can be approximated by 

, ( 2 ) 

a 26A 0 

where SA 0 = [A 2 - Ai] 0 is the local z = 0 separation between both 
wavelengths, and e = 6A z /( 1 + z) - SA Q is the difference between 
the measured line separation at redshift z in rest frame and the local 
one. Thus, in principle, the larger the difference between the pair of 
lines, the better the precision for measuring Aar/a. 

Concerning emission lines, the most suitable pair of lines is 
the [O iii] doublet, which is often present in quasar spectra with rel¬ 
atively high-S/N. The vacuum values for the [O iii] doublet wave¬ 
lengths are 

dj° ‘"I = 4960.295 A 4° m] = 5008.240 A (3) 

<j 4° iii] =47.945 A, (4) 

which are published in the NIST Atomic Spectra Database. 1 These 
transitions are forbidden (they correspond to magnetic dipole and 
electric quadrupole transitions), and they are not observed in the 
laboratory. The wavelength experimental values are obtained indi¬ 
rectly by first computing the energy levels from observed wave¬ 
lengths using a theta-pinch discharge (Pettersson 1982). The wave¬ 
length separation has directly been measured in the infrared from H 
n regions using a balloon-borne telescope and Michelson interfer¬ 
ometer (Moorwood et al. 1980). Both measurements of the wave- 

1 http://physics.nist.gov/PhysRefData/ASD/lines_form.html 


length separation, indirectly with the theta-pinch discharge and di¬ 
rectly with the Michelson interferometer, are in good agreement, 
being the Michelson interferometer more accurate with an error 
< 5 x HT 4 A. 

From equation (2), a determination of e with a precision of 
lA allows for an uncertainty of 10 -2 in A a/a when using the 
[O iii] doublet. The precision from the NIST atomic data allows 
for a determination of A a/a up to 1CT 5 , which is a bit less than 
the uncertainty in our result. One could perform a blind analysis 
in order to search for a possible variation on a, where the absolute 
wavelength values are not required, if one had a large enough sam¬ 
ple distributed in redshift. However, the precision on the absolute 
wavelengths limits the usefulness of high-resolution spectroscopy 
until better measurements of the [O iii] lines (or just their separa¬ 
tion) are available. 


3.2 Implementation 

The code developed for the analysis of the quasar spectra fol¬ 
lows the one described in Gutierrez & Lopez-Corredoira (2010), 
although there are some modifications and more information has 
been extracted from the analysis. We describe the main character¬ 
istics of our code below. 

3.2.1 Wavelength sampling 

We consider only the experimental data together with their errors 
as processed by the SDSS pipeline to obtain the constraint on the 
possible variation of a. We do not resample the wavelength range 
by using an interpolation method. Since the pixel spacing is uni¬ 
form in log-wavelengths, a given range of wavelengths in rest frame 
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Figure 4. Seventh-order polynomial fits (red) to the continuum spectrum with their residuals for three typical quasar spectra at different redshift. The gaps in 
the residuals are the masked regions corresponding to (from right to left) Her, the [O hi] doublet, H/J, H y. 11 6 and Mg n (black dashed lines). 


(/L, /t + ) has the same number of pixels N, i.e. 

J rvt+ (l+z) 3/1+71 

d (log /l) = log -i—-- = log -± , (5) 

•Mt+z) A-(l+z) A- 

and is independent of the redshift of the object. All the wavelength 
intervals with the same width in rest frame will have the same num¬ 
ber of experimental points. 

3.2.2 Fit of the continuum spectrum 

First, we fit a seventh-order polynomial to subtract the continuum 
spectrum while masking regions where strong and wide emission 
lines are present (Hay H/3, Hy, H<5, Mg ii and the [O m] doublet). 
Our method differs from Gutierrez & Lopez-Corredoira (2010) in 
that they use a cubic local spline to fit the continuum masking 
strong emission lines. The chosen order of the polynomial provides 
enough degrees of freedom to reproduce different continuum fea¬ 
tures. In Section 3, we test how our measurement for A ce/a is af¬ 
fected by changing the polynomial order. Hundreds of continuum 
spectra fits were checked by eye. The residuals from the fits are 
smaller than the errors on the flux densities. Fig. 4 shows three dif¬ 
ferent spectra with their continuum fit and residuals. 

3.2.3 Signal-to-noise ratio 

We follow Gutierrez & Lopez-Corredoira (2010) for the determina¬ 
tion of S/N. Hence, we compute the standard deviation of the flux 
between 5040(1 + z) and 5100(1 + z)A (where z is the redshift 
of the quasar) where there are no strong emission or absorption 
lines. Then, we search for the maximum of the [O iii] 5008 line, 
and determine S/Nj-q m j 5008 as the ratio between the maximum of 
the line and the previously computed standard deviation. Although 
for a more reliable determination of the S/N, it is better to use a 
Gaussian fit to the line. This procedure avoids possible issues re¬ 
lated when fitting data with very low S/N. This S/N is used in the 
criterion ii (Section 2) to build our fiducial sample. 

3.2.4 Measurement of the emission-line wavelengths 

To measure the wavelengths of the [O iii] doublet, our fitting code 
needs as input an accurate estimate of the redshift of the quasar, 
at least with an error Az < 3 x 10 -3 . This allows a search for the 


emission lines in a 15 A window around the expected location of 
the [O iii] lines. The SDSS pipeline provides a determination of 
the redshift based on a x 2 fit to different templates; we refer to 
Bolton et al. (2012) for more details. These redshift estimates have 
errors between 1(L 4 and 10 -5 , which are sufficient for our pur¬ 
poses. Moreover, there is also a visual redshift estimation which 
can be found in the quasar catalogue DR12Q (Paris et al. 2015). 
The difference between both redshift estimates (if any) is usually 
kvis - Zpipel ~ 5 x 1CL 4 . We decided to adopt the visual redshifts. 

The centroid positions of the [O iii] emission lines are deter¬ 
mined by four different methods. 

(i) Gaussian profile method. 

First, we search for the maximum flux value in an ~ 15(1 + 
z)A window around the expected position of the line (according 
to the redshift provided by the DR12Q catalogue). This procedure 
automatically erases any bias produced by the redshift value. Then, 
we make an initial Gaussian fit around the position of the maximum 
flux value using a fixed width of ~ 10(1 + z) A. From this first fit, 
we obtain a new position for the line centroid and a Gaussian width. 
These values are used as initial parameters for the final fit of the 
lines; namely, the wavelength range considered to perform the final 
fit is centred around the position of the line centroid, and it is four 
times the Gaussian width of the lines. This approach means that we 
consider pixels up to 2a away from the centre of the line. Hence, 
some lines are fitted using ~ 4-5 pixels, while others with ~ 15-20 
pixels depending on the line width. The fit takes into account the 
flux errors for each pixel, i.e., we use the ivar column found in each 
spectrum as weights for the fit. Our final centroid measurement for 
each considered line corresponds to the centroid of the Gaussian fit 
done in the last step of the adopted procedure. We also derive an 
error for A a/a using the standard errors for the centre position of 
the Gaussians. This is our main method for measuring a. 

(ii) Voigt profile method. 

Following the same procedure than when using a Gaussian pro¬ 
file, we make the fit with a Voigt profile instead of a Gaussian. More 
precisely, we use a pseudo-Voigt profile which is a linear combi¬ 
nation of a Gaussian and a Lorentzian profile. Then, we have one 
more parameter, i.e. the amplitude of the Lorentzian function, while 
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Figure 5. [O in] (left-hand panel) and [Ne hi] (right-hand panel) lines for SDSS-J121417.80+293143.4, at redshift z = 0.063. The measured A a/a for this 
quasar is Aa/ag u “ j = (2.3 ± 7.6) x 1(T 4 , Aa/a^j^ = (3.3 ± 12.6) x 1(T 4 and Aa/aj^f, = (39 ± 8) x 1(T 4 , Aa/aj^' = (37 ± 9) X 1(T 4 . The 
measured Aa/a for [Ne hi] is not consistent with zero regardless of the profile; see Fig. 14 and last paragraph of Section 4 for discussion. Each panel shows 
the flux density for each pixel with their respective error bars (solid symbols), together with the Gaussian fit (dotted red curve) and the pseudo-Voigt profile 
(thick grey curve) to each of the lines. The fitting procedure (described in the text) only takes into account the experimental data (solid symbols) weighted by 
their error bars. Notice how the deviation of the line centroid position derived from our Gaussian fit (vertical solid line) with respect to the expected position 
of the line (vertical dashed line) according to the visual redshift provided by the DR12Q catalogue are well correlated for the same pair of [O hi] and [Ne iii] 
lines, and for the different set of lines. The green shaded vertical areas highlight the uncertainty for the expected position of the lines due to the quasar redshift 
error (» 5 X 10 4 ). Also shown is a fourth-order spline interpolation to the spectrum after subtracting the continuum (thin solid line). The [Ne iii] lines are 
weaker by one order of magnitude than the [O iii] lines, which is usually the case for all the spectra showing both pair of lines. The weak line near the stronger 
[Ne iu] line is blended with He i (3889.75 A) and Htf (3890.16 A). 


its width and its position are the same as those for the Gaussian 
profile. 

In Fig. 5, we depict the [O iii] and [Ne iii] lines for the 
same quasar spectrum to illustrate the Gaussian and Voigt fitting 
methods. 

(iii) Integration method. 

Here, the centroids of the lines are obtained by integrating 
around l<x from the position of the fitted Gaussian profile. This 
technique provides indications of whether there is H/3 contamina¬ 
tion. However, due to the mid-resolution of the spectra R * 2000, 
this method is not very accurate. 

(iv) Modified Bahcall method. 

In Bahcall et al. (2004) the authors used a different approach 
to compute the line positions. They performed a third-order spline 
interpolation to the stronger [O iii] 5008 line, then fitted this inter¬ 
polation to the weaker 4960 line by adjusting the amplitude and 
separation of the profile. We have modified this method by using a 
Gaussian fit to the stronger line rather than a third-order spline. 

Although we have described four different methods, the main 
results for Aa/a presented in this work are based on the Gaussian 
fitting method, while the other three are used only for comparison 
(see Section 4). 

Finally, our final result for Aa/a and its error is obtained in the 
same way as in Chand et al. (2005), namely we compute a weighted 
mean and a weighted standard deviation, where the errors for Aa/a 
of each spectrum are used as weights. 


3.3 Simulated spectra 

In order to test the robustness and accuracy of our methodology, 
we generate realizations of quasar spectra using as noise a nor¬ 
mal distribution centred at the flux value, and taking the error in 
each pixel as the standard deviation. From our fiducial sample 
(10 363 quasars), we simulate 100 realizations for each spectrum 
(> a million in total). This number of realizations provides rea¬ 
sonable statistics to derive an error from the standard deviation of 
the measurements on the realizations of each real spectrum, while 
the computation time remains reasonable (~ 2 d) using a standard- 
size computer. The estimated error derived from the simulations 
A(Aa/ar) s i m includes 

A(A a/a) 2 sim = A(Aa/af 6l + A(Aa/a)j ontinuum + A(A a/af coAt , (6) 

where A(Aa-/o') flt is the error derived from the Gaussian fits, which 
is our error estimate for each real spectrum; AfAcr/oOcontinuum is 
the error from different continuum subtraction due to the Gaussian 
noise, and A(A(r/tt) code is the systematic error of our code. Then, 
we expect AfAcr/alsim > A(AQ-/a)j t and their difference will be an 
indication of the continuum and systematic errors. 

Fig. 6 (left-hand panel) shows the correlation between the er¬ 
ror in Aa/a from the Gaussian fits of each real spectrum and the 
standard deviation for Aa/a of its 100 realizations. The standard 
deviations from the simulations are within a factor of 0.5 - 2 of the 
standard errors from the fits for 97% (84%) of the cases when both 
quantities are < 5 x 10~ 3 (< 50 x 10~ 3 ). This shows that our code 
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A(Aof/or)(i t (x 10 -3 ) Redshift 

Figure 6. Left-hand panel: errors for Aa/a obtained from the simulations (standard deviation of the Aa/a measurements on 100 realizations of each real 
spectrum) and standard errors from the Gaussian fits for our fiducial sample. The solid line represents a one-to-one correspondence, while the dashed lines 
have slopes of 2 and 0.5. Only the simulation and fit errors smaller than < 5 X 1 0 3 are shown. Right-hand panel: errors estimated from the simulations as a 
function of redshift. Spectra with A ( Aa/a ) flt > 5 X 1 0 1 are shown as red crosses (24% of the total). There is a clear division between two different set of 
spectra which correlates with the sky emission lines (see discussion in the main text). 



-10 -S 0 5 10 0 2 4 6 8 10 


(Atr/tl) Qaussian (8 10 ) Error (Ao'/o') Gaussian (8 10 ) 

Figure 7. Left-hand panel: measurements of Aa/a using Gaussian and Voigt fitting profiles. Non-compatible measurements at lcr are shown as red crosses 
(6.5% of the total). Right-hand panel: errors from the Gaussian and Voigt fitting. Non-compatible measurements at I <r are shown as red crosses. 


and the continuum subtraction do not introduce noticeable system¬ 
atic errors compared to the Gaussian fitting. However, there is a set 
of spectra (9% of the total) for which the simulations provide much 
larger errors A(A a/a) >0.1. Fig. 6 (right-hand panel) shows the 
errors from the simulations as a function of redshift for our fidu¬ 
cial sample. Red crosses stand for spectra whose Gaussian fit error 
A(Aa/ar) fl , > 5 x 10~ 3 (24%). The errors are distributed in two 
clouds of points. For the cloud with A(Aa/af) s i m ~ 1, the virtual 
realizations of each spectrum seem to differ significantly from the 
real spectrum. Since we use the error in each pixel to build the re¬ 
alizations, the relative error is large for these spectra, which is an 
indication of a low S/N ratio or large absolute errors in the pixels, 
for instance in wavelength regions with sky emission lines. In fact, 
the cloud with bigger errors mimics the sky spectrum. Note also the 
strong correlation between this cloud of points and the spectra with 
large Gaussian fitting errors (red crosses). The other set of points 
with A(Ao7a) sim ~ 10 -3 are close to our error estimation on the 
measurement of Aa/a based on the Gaussian fits. 


As a further proof, we also simulate realizations of the 1416 
dropped spectra because of sky emission lines (criterion iv, see 
Section 2). We found that more than 80% of the spectra have 
AfAtt/oOsim > 0.1. This confirms that these spectra have very low 
S/N and/or large pixels error due to the proximity of the lines to 
strong sky emission lines. 

3.4 Gaussian versus Voigt fitting profiles 

The results obtained when using Voigt profiles instead of Gaussian 
ones are compared in Fig. 7. The Voigt and Gaussian measurements 
are lcr-compatible for the 93.5% of the cases (98.3% at 2<x). Re¬ 
garding the errors, there is no clear improvement when using either 
of both methods. However, Voigt profiles have one more parameter 
and restrict the number of degrees of freedom. Due to the spectral 
mid-resolution and the fact that the [O m] lines are very narrow, 
there are often only a few pixels to fit, which frequently lead to non- 
convergent fits. This reduces the quasar sample in « 1000 quasars. 
Further discussion about both profiles can be found in Section 4. 
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4 SYSTEMATICS 

In this section, we examine the possible unnoticed systematic er¬ 
rors by analysing different quasar samples. Table 2 summarizes all 
the samples considered together with their mean redshifts and the 
measured value for A a/a. 

We consider the following sources of systematic errors. 

(i) Misidentification of the lines. The expected line widths and 
amplitudes are useful to avoid misidentification of the [O in] emis¬ 
sion lines, (a) Line widths: since both lines originate on the same 
upper energy level, their width must coincide. We check that this 
is the case by considering quasars whose [O hi] line widths are the 
same within a relative fraction. For more than half of our fiducial 
sample, the [O m] line widths differ by less than 10% (see Table 2). 
(b) Amplitude ratio: atomic physics states that the amplitude ratio 
between the [O m] 5008 and [O m] 4960 lines is 2.98 (Storey & 
Zeippen 2000) (as quoted in Section 5, we obtain 2.96 ± 0.02 syst ). 
Thus, we consider different samples where this ratio differs by less 
than a certain amount from 2.98 (see Table 2). All the samples con¬ 
sidered in this test yield results for A a/a compatible with zero. 
Fig. 8 displays the Gaussian widths and fluxes of both [O in] emis¬ 
sion lines for our fiducial sample. 

(ii) Windows for the Gaussian fits. We use a wavelength range 
of 2cr around each [O m] line in order to obtain the final Gaus¬ 
sian fit to the line profiles. We study how our results depend on 
this choice. By considering a larger wavelength interval, the results 
are more affected by the H/f contamination and possible asymme¬ 
tries on the line wings. The differences in the number of spectra for 
these samples [which are obtained by applying the selection crite¬ 
ria (i)-(iv) discussed in Section 2.1] arise because of the criteria 
concerning the non-converging fits and the sky emission lines de¬ 
scribed in Section 2. 

(iii) H/l contamination. We analyse samples where the ratio be¬ 
tween S/Nh/j and S/Nj-q m j 4%0 is constrained. Despite the fact that 
the value for Aa/a decreases as we place more stringent constraints 
on Hyfl, it is always consistent with no variation in a within the er¬ 
rors. This analysis demonstrates that the strength and/or width of 
the H/3 line do not affect substantially the result for A a/a when a 
weighted mean is adopted. 

(iv) Continuum subtraction. We use a seventh-order polynomial 
to subtract the continuum spectrum. We examine if the polynomial 
order has important effects on our measurements. Our values for 
A a/a and their errors are only slightly affected by the chosen poly¬ 
nomial order. 

(v) Goodness of Gaussian fits. We quantify the quality of the 
Gaussian fits by the R 1 coefficient. All the considered samples show 
values for A a/a consistent with no variation in a. 

(vi) Broad lines. We also study samples where the width of both 
lines is less than a certain value (in kms 4 ). These samples are 
consistent with no variation of a. Samples built from narrow lines 
< 300 kms -1 may be more affected by misidentification of [O hi] 
lines as sky lines. 

(vii) Different methods for measuring the [O m] line position. 
We compare the results obtained by the methods to measure the 
position of the [O iii] lines described in Section 3.2.4. Since not 
all the methods provide an error for the measurement, we cannot 
calculate a weighted mean, and it is necessary to select a more re¬ 
stricted sample. Then, we consider a sample where the difference 
between the widths of the lines is less than 25%, the amplitude ra¬ 
tio is constrained to differ from the theoretical value 2.98 (Storey & 
Zeippen 2000) by less than 0.5, and the S/N H /j is smaller than half 
the S/N [0 m] 4960 - 


Table 2. Results for Aa/a considering several samples with different con¬ 
straints. The number of quasar spectra, the mean and standard deviation of 
the redshift and the value for Aa/a are shown. 


a 4960/5008 - 1 

No. of quasar spectra 

Redshift 

A a/a (xlO -5 ) 

< 50% 

10028 

0.56 ±0.21 

1.6 ±2.3 

< 25% 

8877 

0.56 ±0.21 

1.9 ±2.3 

< 10% 

5846 

0.56 ±0.21 

1.7 ±2.5 

< 5% 

3458 

0.54 ±0.22 

-0.9 ± 3.0 

[FA x Cr] 5 008/4960 

No. of quasar spectra 

Redshift 

Aa/a (xlO -5 ) 

2.98 ± 0.50 

8327 

0.56 ±0.21 

1.8 ±2.4 

2.98 ± 0.25 

5761 

0.55 ±0.21 

-0.4 ± 2.6 

2.98 ±0.10 

2658 

0.54 ±0.21 

0.0 ± 3.4 

2.98 ± 0.05 

1411 

0.52 ±0.22 

5.2 ±4.6 

Fit width 

No. of quasar spectra 

Redshift 

Aa/a (xlO -5 ) 

2 cr 

10 363 

0.56 ±0.21 

1.4 ±2.3 

3 cr 

10252 

0.59 ± 0.20 

5.5 ±2.5 

4cr 

9978 

0.59 ±0.20 

7.1 ±2.7 

5cr 

9726 

0.59 ±0.20 

5.3 ±2.6 

S/Nn^/[Oiii]4960 

No. of quasar spectra 

Redshift 

A a/a (XlO -5 ) 

< 5 

10 338 

0.57 ±0.21 

1.4 ±2.3 

< 2 

9831 

0.57 ±0.21 

0.6 ±2.3 

< 1 

8162 

0.57 ±0.21 

0.1 ±2.5 

< 0.5 

5831 

0.58 ±0.21 

-0.7 ±2.8 

Pol. order (cont.) 

No. of quasar spectra 

Redshift 

Aa/a (XlO -5 ) 

3 

10528 

0.57 ±0.21 

1.0 ±2.3 

5 

10550 

0.57 ±0.21 

1.3 ±2.3 

7 

10 363 

0.56 ±0.21 

1.4 ±2.3 

9 

10471 

0.56 ±0.21 

-1.1 ±2.3 

R 2 (both fits) 

No. of quasar spectra 

Redshift 

Aa/a (XlO -5 ) 

> 0.9 

9254 

0.56 ±0.21 

1.5 ±2.4 

> 0.97 

6045 

0.56 ±0.21 

2.8 ±2.7 

> 0.99 

2301 

0.54 ± 0.21 

2.0 ±3.5 

> 0.995 

845 

0.51 ±0.22 

-0.4 ± 4.8 

[O iii] 5008 (kins No. of quasar spectra 

Redshift 

A a/a (xlO -5 ) 

< 1000 

10 353 

0.56 ±0.21 

1.4 ±2.3 

< 500 

8990 

0.56 ±0.21 

0.2 ±2.4 

< 300 

2798 

0.52 ±0.22 

-6.8 ±3.9 

< 200 

150 

0.52 ±0.24 

21 ± 18 

Method 

No. of quasar spectra 

Redshift 

Aa/a (XlO -5 ) 

Gaussian (weighted) 

4537 

0.58 ±0.20 

-0.4 ± 2.8 

Gaussian 

4537 

0.58 ±0.20 

1.2 ±4.5 

Integration 

4537 

0.58 ± 0.20 

3.6 ±4.8 

Modified Bahcall 

4537 

0.58 ±0.20 

0.8 ±4.4 

Median 

4537 

0.58 ±0.20 

1.8 ± 1.4 

Gauss versus Voigt 

No. of quasar spectra 

Redshift 

Aa/a (XlO -5 ) 

Gaussian profiles 

8485 

0.55 ±0.19 

0.4 ± 2.5 

Voigt profiles 

8485 

0.55 ±0.19 

-1.1 ±2.8 

Mixed profiles 

8485 

0.55 ±0.19 

1.3 ±2.4 


(viii) Gaussian versus Voigt profiles. We compare the results 
for 8485 quasars from our fiducial sample after dropping 1878 
spectra with non-converging Voigt fits (this reduction increases 
the statistical error). We also compute a ‘mixed’ value for Aa/a 
where for each spectrum we use the value for the variation of 
the fine-structure constant with smaller error, either (Aff/a) Gauss or 
(A a / O') Voigt- 
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Figure 8. Left-hand panel: Gaussian widths (in kms -1 ) for both [O in] lines. Both lines originate on the same upper energy level, then their widths must 
coincide (red dashed line). Right-hand panel: fluxes for both [O in] lines. The theoretical flux ratio is 2.98 (red dashed line). The entire fiducial quasar sample 
is shown. 
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Table 3. Detailed information about the bins in Fig. 10. 
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Figure 9. Errors for A a/a derived from the Gaussian fits (grey points) for 
our fiducial sample, moving mean of the these errors (blue line) using over¬ 
lapping bins (100 spectra per bin, Az ~ 0.025), moving standard deviation 
of A a/a measurements using the same bins (red line) and a typical sky 
spectrum. 


We have also analysed the standard deviation and errors of 
the results for A a/a as a function of redshift (Fig. 9). Even though 
we have imposed a constraint on our initial sample based on the 
sky emission lines, the standard deviation and errors still correlate 
with the sky. In particular, for the correlation with the moving stan¬ 
dard deviation, this means that the precision in our measurement of 
Aa/a along the whole redshift interval is limited by the sky sub¬ 
traction algorithm. 

5 RESULTS 
5.1 [O m] lines 

We used a total of 10 363 quasar spectra, drawn from the SDSS- 
III/BOSS DR12Q catalogue, after applying the selection criteria 
(i)-(iv) (see Section 2), to measure the possible variation of the 
fine-structure constant. The following measurement is obtained: 

— = (1.4 ± 2.3) X 10 -5 . 
a 


Redshift interval 

No. of quasar spectra 

Redshift 

A a/a (xl0~ 5 ) 

0.390 - 0.460 

817 

0.42 ± 0.02 

-5.2 ±6.8 

0.460 - 0.520 

723 

0.49 ± 0.02 

5.5 ± 8.9 

0.520 - 0.580 

757 

0.55 ± 0.02 

0.4 ± 9.2 

0.580 - 0.625 

843 

0.60 ± 0.01 

40.4 ± 9.4 

0.625 - 0.675 

988 

0.65 ± 0.01 

-3.5 ± 7.4 

0.675-0.715 

1299 

0.69 ± 0.01 

-8.2 ±7.4 

0.715-0.765 

1117 

0.74 ± 0.01 

1.7 ± 7.1 

0.765 - 0.820 

1444 

0.79 ± 0.02 

18.1 ± 8.3 

0.820 - 0.880 

644 

0.84 ± 0.02 

4.7 ± 9.2 

0.880- 1.000 

580 

0.93 ± 0.03 

17.0 ± 13.3 


This value is consistent with the previous results reported in differ¬ 
ent investigations based on the same method: Bahcall et al. (2004), 
Gutierrez & Lopez-Corredoira (2010), and Rahmani et al. (2014). 
The redshift dependence of the measurements is shown in Fig. 10 
(left-hand panel), where several bins have been made taking into 
account the redshift intervals affected by the sky (shaded zones). In 
the right-hand panel, we show the results obtained from the simu¬ 
lations described in Section 3, using the same redshifts intervals for 
the bins. The main differences between the real results and the sim¬ 
ulations are in the regions where there are strong sky lines (shaded 
regions), while being in agreement in the remaining zones. Detailed 
information about each bin for the real data can be found in Table 3. 

Our results are little affected by the specific constraints im¬ 
posed in our sample as discussed in Section 4. For instance, we 
vary the width for the Gaussian fits, the contamination of H fi, the 
polynomial order used to fit the continuum spectrum, the quality of 
the Gaussian fits and test different methods to measure A a/a. The 
most important effect found is that by considering broader widths 
for the Gaussian fits, the results are more affected by the contamina¬ 
tion from HyS and possible asymmetries in the line wings. We have 
also checked for possible misidentifications of the [O hi] emission 
lines using their expected widths and amplitude ratio. 

Table 4 contains the results for A a/a when the lower bound 
on the S/N[Q m j 500g is increased. All the results remain consis¬ 
tent with no variation of the fine-structure constant. In Fig. 11, 
the measured Aa/a for our fiducial sample as a function of the 
S/N[ 0 m ] 5008 are plotted together with their errors. 
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Figure 10. Left-hand panel: A a/a versus redshift (real data). Details about each bin are listed in Table 3. Right-hand panel: (Aa/a) sim versus redshift 
(simulations). A typical sky spectrum and shadowed regions where the sky contamination is expected to be large, are shown as reference. 


Table 4. Results for several samples with different constraints on the 
S/N[o hi] 5008 ■ I ' 01 eac h sample, the number of quasar spectra, the mean 
redshift, together with its standard deviation and the value for Aa/a are 
shown. 


S/N[0 III] 5008 

No. of quasar spectra 

Redshift 

Aa/a (xlO -5 ) 

> 10 

10 363 

0.56 ±0.21 

1.4 ±2.3 

> 20 

5270 

0.53 ± 0.21 

-0.5 ± 2.5 

> 50 

1498 

0.47 ± 0.20 

-3.4 ±3.1 

> 100 

451 

0.41 ± 0.19 

-2.0 ± 3.6 

> 500 

12 

0.24 ±0.19 

6 ± 12 


Table 5. Results for the North and South galactic hemispheres. 


Galactic hemisphere No. of quasar spectra 

Redshift 

Aa/a (xl0~ 5 ) 

North 

8069 

0.56 ±0.21 

2.6 ±2.6 

South 

2294 

0.59 ± 0.20 

-3.1 ±4.9 


The distribution of BOSS quasars in the sky (see Fig. 1, left- 
hand panel) suggests to divide the sample into two, one for the 
North galactic cap and one for the South galactic cap. Table 5 de¬ 
scribes the results for each galactic hemisphere, and no statistical 
meaningful difference is found. In order to look for a spatial varia¬ 
tion, we also carried out a more precise analysis by fitting a dipole. 
First, we fixed the direction on the sky of the dipole and performed 
a linear fit (A a/a(8) = a cos 6 + b) of the measurements of the vari¬ 
ation of the fine-structure constant as a function of cos 6, where 9 is 
the angle between the dipole axis and a quasar in the sky. Different 
fits were done for the dipole axis lying in a grid in RA and Dec. 
(~ 1° x 1°). However, there is not statistical significance to deter¬ 
mine the dipole axis with a meaningful error, i.e. smaller than the 
whole sky. There has been a claim on a significant deviation of a 
from being a constant as a function of space (King et al. 2012), with 
a dipole amplitude 0.97^-^ xlO -5 in the direction RA= 17.3±1.0h 
and Dec.= -61° ± 10°. Fixing the dipole in that direction, we get 
(-4.3 ± 3.4) x 1CT 5 for the dipole amplitude and (0.8 ± 2.3) x 1(T 5 
for the monopole term, which are not precise enough to compare 
with that work. 

We are inclined to parametrize the possible time variation of 
a with redshift z. This is justified since any possible variation on 
a must be dominated by the local geometry of space-time (at least 


20 50 too 200 500 1000 



20 50 too 200 500 1000 


S/A'ioiiqsoos 

Figure 11. Aa/a versus S/Njq n - 0(lx (top panel) with the moving stan¬ 
dard deviation (black lines) using overlapping bins (100 spectra per bin) 
and the error on Aa/a (bottom panel) in linear-log scale for our fiducial 
sample. The deviation of A a/a from zero and its error steadily decreases as 
the S/N increases. 

if we consider the dynamics of the Universe as the main reason 
for such variation). Therefore, one is led to consider the possible 
variation of or as a function of redshift (z = l/a(t) - 1) or the Ricci 
scalar (R(t) = 6H(t) 2 \ I - g(f)J), where a(t) is the scale factor, H(t) 
the Hubble parameter and q(t) is the deceleration parameter. Since 
the Ricci scalar is not known for each quasar, it is straightforward 
to consider a possible variation with redshift. In contrast, for a time 
parametrized model of the variation of a the analysis depends on 
the particular cosmology considered. Since there is no significant 
clear dependence, we use a linear model in redshift. Then, for 

A a/a = az + b , (7) 

we obtain 

a = (0.7 ± 2.1) x 1(T 4 ; b = (0.7 ± 1.4) x 1(T 4 ; (8) 

which do not show any dependence of A a/a with redshift. 
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Table 6. Results for the line ratio when polynomials of different orders 
are used to subtract the continuum and different range for the Gaus¬ 
sian fits are used. For each sample, the number of quasar spectra, the 
mean redshift together with its standard deviation and the value for 
[F,! X cr] 5008 / [F ,1 X c ] 4 %o are shown. 


Polynomial order No. of quasar spectra Redshift [F,i X cr] 5 oo 8/4960 


3 

10528 

0.57 ±0.21 

2.96 

5 

10550 

0.57 ±0.21 

2.94 

7 

10 363 

0.56 ±0.21 

2.96 

9 

10471 

0.56 ±0.21 

2.98 

Fit width 

No. of quasar spectra 

Redshift 

[F.l X (±15008/4960 

2 o- 

10 363 

0.56 ±0.21 

2.96 

3 cr 

10252 

0.59 ± 0.20 

2.92 

4 cr 

9978 

0.59 ± 0.20 

2.98 

5cr 

9726 

0.56 ±0.21 

2.98 


Table 7. Results for SDSS-II/DR7 and BOSS (SDSS-III/DR12) samples. 
For each sample, the number of quasar spectra, the mean redshift, together 
with its standard deviation and value for A a/a are shown. 

Sample 

No. of quasar spectra 

Redshift 

Aa/a (xl0~ 5 ) 

DR7 

2853 

0.38 ±0.15 

0.5 ± 2.8 

DR7 (SDSS cont.) 

3009 

0.38 ±0.15 

-0.4 ± 2.7 

BOSS (DR12) 

10 363 

0.56 ± 0.21 

1.4 ±2.3 

BOSS + DR7 

13175 

0.51 ± 0.21 

0.9 ± 1.8 


From this sample, we also obtain a value for the line ratio 
[Fj x tr] 5 oo 8 / [F A x cr] 4 96 o = 2.96 ± 0.02 syst , where F A is the max¬ 
imum flux density of the line, and a is the Gaussian width. The 
value reported is a weighted mean where the S/Nj-q m j 500g is used 
as weights. The quoted systematic error is computed from the anal¬ 
ysis of samples with different polynomial orders for the continuum 
fit and different widths for the line fitting (see Table 6), since this 
quantity is more affected by these two parameters. The value we 
obtain is in agreement with the best current theoretical value, i.e., 
2.98 (Storey & Zeippen 2000). 

Finally, we have also considered quasar spectra from the 
SDSS-II/DR7, which were observed using the previous spectro¬ 
graph instead of the upgraded BOSS spectrograph (see Section 2). 
From the DR7 quasar catalogue (Schneider et al. 2010), which 
contains 105 783 quasars, we select a sample of 2853 quasars 
up to redshift z = 0.8 using the same criteria described in Sec¬ 
tion 2. This number is similar to the quasar spectra considered 


Table 8. Values of Aa/a using the combined sample BOSS+DR7 for the 
same redshift bins as in Table 3. 


Redshift interval 

No. of quasar spectra 

Redshift 

Aa/a (xl0~ 5 ) 

0.390 - 0.460 

1279 

0.42 ± 0.02 

-2.5 ± 4.8 

0.460 - 0.520 

1076 

0.49 ± 0.02 

7.2 ± 6.4 

0.520 - 0.580 

1071 

0.55 ± 0.02 

1.1 ±7.1 

0.580 - 0.625 

1025 

0.60 ± 0.01 

30.8 ± 8.1 

0.625 - 0.675 

1191 

0.65 ± 0.01 

-5.1 ±6.5 

0.675 - 0.715 

1424 

0.69 ± 0.01 

-6.2 ± 7.0 

0.715-0.765 

1220 

0.74 ± 0.01 

1.4 ±6.8 

0.765 - 0.820 

1519 

0.79 ± 0.02 

15.0 ±8.1 

0.820 - 0.880 

644 

0.84 ± 0.02 

4.7 ± 9.2 

0.880- 1.000 

580 

0.93 ± 0.03 

17.0 ± 13.3 



Figure 12. Difference between (Aq , /q')boss and (Aa/a/ota measurements 
for the 41 quasars observed by SDSS-I/II that were re-observed by BOSS. 
Both values for Aa/a are consistent. 

by Rahmani et al. (2014, Table 1). We re-analyse this sample us¬ 
ing the methodology presented in this work, and we find Aa/a = 
(0.5 ± 2.8) x 10~ 5 . By combining this DR7 sample with our fiducial 
BOSS (DR12) quasar sample (after eliminating 41 spectra which 
were re-observed by BOSS), we obtain our final sample which con¬ 
tains a total of 13 175 quasars. The value obtained for this combined 
sample is reported as a final result of this investigation: 

— = (0.9 ± 1.8) x 1(T 5 . 
a 

Table 7 shows the results for DR7, DR7 using the continuum fit 
provided by the SDSS pipeline, 2 BOSS (DR 12) and the combined 
BOSS+DR7. It can be seen that the mean redshift for the DR7 sam¬ 
ple is lower than that for BOSS. Note that there is also a big differ¬ 
ence on the mean S/Nj-q m j 500g of both samples: S/Nj^ 7 ^ = 

60 and S/N^S 88 . = 33, which also explains why the statisti- 

[O III] 5008 r 

cal errors for Aa/a do not reflect the expected reduction due to the 
increase in the number of quasars of our BOSS sample. Table 8 
shows the results of Aa/a using the combined sample in the same 
redshift bins considered for our fiducial sample. Fig. 12 shows the 
difference of the values obtained for A a/a for the 41 re-observed 
quasars. Both BOSS and DR7 measurements are in perfect agree¬ 
ment within the error bars. This test is a good check for the relia¬ 
bility of our code and the consistency of the SDSS spectra obtained 
with different spectrographs. 

There are massive galaxy surveys to be carried out during the 
next decade. For instance, eBOSS and DESI will take spectra from 
millions of galaxies. Therefore, it is interesting to give an estima¬ 
tion of the accuracy that will be obtained when using galaxy spectra 
instead of quasars. For this, we have analysed the galaxy spectra 
collected by the DEEP2 survey (Newman et al. 2013) taken with 
resolving power = 6000. From this sample, we found 4056 galax¬ 
ies with strong [OIII] lines. Naively, one would expect that the error 
on Aa/a should be 

A(AQ'/a) galaxics DEE p2 ~ /sample X finst X /object X A (AQ'/Q') quasai . s BOSS , 

where / sa mpie(= V10 363/4056) takes into account the different 

2 The SDSS pipeline provides a continuum fit for the DR7 spectra. The 
good agreement between the value for Aa/a obtained with the SDSS con¬ 
tinuum fit and our continuum fit is a good test for our code. 
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Figure 13. Errors on A a/a as a function of the FWHM of [O in] 5008 for 
our fiducial BOSS quasar sample (blue points) and for the DEEP2 galaxy 
sample (red crosses). 

number of objects in each sample, / inst (=s 2000/6000) stands for 
the different resolution of the spectra and /object is an extra factor 
due to the different characteristics of quasar and galaxy emission 
lines which affect the uncertainty of the line positions. This last 
factor is proportional to the line widths and inversely proportional 
to the line fluxes: 

r FWHM ga , axies ^ / FIuXg a | aMes j gins 

/objeCt “ FWHM qllasars X \ Flux quasars / ' 

These numbers for the [O in] 5008 line are approximately 
FWHMgaiaxies at 120 kms" 1 , FWHMq Uasars at 420 kms -1 , 
Flux galaxies at 70 and Flux quasars at 210 in units of 10 -17 erg cm -2 s -1 , 
obtained from the DEEP2 sample and from our fiducial sample. 
Thus, the expected error is 1.1 x 10 5 , where we have considered the 
error of our fiducial sample A (Aa/ar) quasarsBOSS = 2.3 x 10 -5 . Ap¬ 
plying the same criteria described in Section 2 to the DEEP2 galaxy 
sample, we get Aa/a = (-0.9 ± 1.6) x 10 -5 . Thus, the upcoming 
future galaxy surveys will be quite competitive for constraining the 
variation of the fine-structure constant at low redshift z <2. Fig. 13 
shows the error on Aa/a for the DEEP2 and BOSS samples. 

5.2 [Ne m] lines 

We also measure from 462 quasar spectra with [Ne m] emission 
lines the following constraint on the fine-structure constant: 

Aa/ff[Ne In ] = (34 ± 1) X 10 -4 , (11) 

to be compared with 

Aa/ff[Nem] = (36±l)xl<r 4 (12) 

obtained by Gutierrez & Lopez-Corredoira (2010). The analysis 
of the [Ne hi] lines reveals the same systematic effect previously 
observed, namely a clear tendency for a positive variation of a. 
Fig. 14 compares the results obtained for Aa/a for spectra where 
both [O hi] and [Ne hi] lines are present. To account for this effect, 
a shift ~ 0.6 A on the theoretical or observed values of the wave¬ 
lengths for the [Ne iii] lines is necessary. There are experimental 
(Bowen 1955) and indirect (Kramida & Nave 2006) values for the 
wavelengths of the [Ne iii] lines which are in agreement with errors 
at 3 x 10 -2 A. We use the NIST values for the [Ne iii] lines 

2[Ne m ] = 3869.86 A 4 Ne IIl] = 3968.59 A (13) 


Figure 14. [Ne iii] and [O iii] measurements for A a/a. Empty symbols 
stand for spectra with S/N^ e In j 3869 < 35 and solid squares for spectra 
with S/N[N e m ] 3869 > 35. [Nem] measurements have a clear tendency 
to a positive variation of a, which is due to a systematic effect affect¬ 
ing the [Ne iii] measurement. The same effect has already been noticed by 
Gutierrez & Lopez-Corredoira (2010), and it is explained in Section 5.2. 

SA^ = 98.73 A. (14) 

The results for the [O in] doublet guarantee the good calibration 
of the SDSS spectra (and many more independent scientific re¬ 
sults based on the SDSS spectra). Thus, we have measured the 
[Ne in] lines using a high-resolution optical spectrum from the 
planetary nebula IC 418. The IC 418 optical spectrum (3600-7200 
A) was taken under service time at the Nordic Optical Telescope 
(NOT; Roque de los Muchachos, La Palma) in 2013 March with 
the FIES spectrograph. We used FIES in the low-resolution mode 
(R at 25 000) with the 2.5 arsec fibre (centred at the central star 
of IC 418). Three exposures of 1200 s each were combined into a 
final IC 418 spectrum, reaching a S/N (in the stellar continuum) of 
~60 at 4000 A and in excess of ~ 150 at wavelengths longer than 
5000 A (see Dfaz-Luis et al. 2015, for more observational details). 
To measure Aa/a, we need to know the ratio 

K = [W 2 -*)/W 2 + *)] 0 , (15) 

which is independent of the peculiar velocity of the planetary neb¬ 
ula. From our data, we obtain 

R = (1259561 ±4)x 10 -8 , (16) 

compared to the one using NIST values for the wavelengths 

Knist = 1259560 x 10 -8 . (17) 

The difference between the two values translates into a variation 
on Aa/a < 10 -6 . Thus, the measured wavelength separation for 
the [Ne m] doublet does not account for the positive variation on 
a observed using these lines. Fig. 15 (left-hand panel) shows the 
Gaussian fit to the [Ne iii] line profiles present in the IC 418 spec¬ 
trum. 

The IC 418 spectrum shows two different lines near the [Ne hi] 
3968 A line (see Fig. 15, right-hand panel). The stronger one is 
He 3971 A, the other one is He i 3965 A. Hence, we search for 
a possible blending of the [Ne m] line 3968 with these two lines 
in our much lower spectral resolution quasar spectra. Fig. 16 (left- 
hand panel) shows stack quasar spectra with broad [Ne iii] emission 
lines. It can be seen that the weak [Ne iii] line is blended. 

To quantify the displacement produced by the blending with 
He line, we did a Gaussian convolution of the Planetary Nebula 
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Figure 15. Left-hand panel: [Ne hi] 3869 and [Ne in] 3968 together with our Gaussian fits (solid line) from a high-resolution (R » 25 000) spectrum of the 
Planetary Nebula IC 418 obtained with the FIES spectrograph at the NOT telescope. Right-hand panel: IC 418 spectrum centred at [Ne in] 3968 line. The two 
close lines are H e 3971 A and He i 3965 A. 
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Figure 16. Left-hand panel: median-stacked quasar spectra with broad [Ne hi] lines (increasing line width from bottom to top). Both [Ne iii] lines are shown 
(AA 3869,3968 A). The weak [Ne in] line is blended with the two lines H e 3971 A and He i 3965 A. Right-hand panel: shift produced by the H e line in the 
[Ne in] 3968 line as a function of the line intensity ratio of both lines as measured from the Planetary Nebula convolved spectrum. 


spectrum to lower the resolution down to R ^ 2000. Since the line 
intensity ratio of [Ne iii] and H e may differ in the quasar narrow 
emission-line region and the Planetary Nebula, we show in Fig. 16 
(right-hand panel) the shift produced by the H e line as a function 
of the ratio [Ne m] /H e. We get a shift ~0.6 A when H e/[Ne m] is 
~0.5. This explains the systematic found when using [Ne iii] lines 
to measure the variation of the fine-structure constant A a/a in pre¬ 
vious studies (Gutierrez & Lopez-Corredoira 2010; Grupe et al. 
2005). 

6 SUMMARY 

The main conclusions of this work are as follows. 

(i) From 45 802 objects at z < 1 classified as quasars in the 
SDSS-III/BOSS DR12 quasar catalogue, we have extracted a sam¬ 
ple of 10 363 quasars with [O iii] emission lines. Combining this 
fiducial sample with a sample of 2853 previously studied SDSS- 
II/DR7 quasars, we got a final sample of 13 175 after eliminating 
41 re-observed quasars. 


(ii) With this combined sample, we have estimated a value 
for the possible variation of the fine-structure constant of 
Aa/a = (0.9 ± 1.8) x 10~ 5 , which represents the most accurate re¬ 
sult obtained with this methodology. 

(iii) We have also studied how much our results change when 
analysing the fiducial sample according to different properties 
(width, amplitude, S/N and R 2 coefficient of the [O iii] lines), and 
when modifying some parameters of the analysis (polynomial order 
for the continuum subtraction, different methods to determine the 
line position, e.g. Gaussian/Voigt profiles). We conclude that our 
results are quite robust, and they are consistent with no variation of 
the fine-structure constant. 

(iv) From over one million simulated realizations of quasar 
spectra, we conclude that the precision of our emission-line method 
is dominated by the error from the Gaussian fits. Flence, the error 
from the continuum subtraction and any possible systematics from 
our code are small. 

(v) The standard deviation of the results as a function of redshift 
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correlates with the sky. This result suggests that our main source of 
uncertainty is determined by the sky subtraction algorithm. 

(vi) We have determined the ratio of the [O m] transition lines 
to be 2.96 ± 0.02 syst , which is in good agreement with previous 
experimental and theoretical values. 

(vii) The same systematic effect previously noticed by Gutierrez 
& Lopez-Corredoira (2010) has been found on the [Ne in] lines 
measurement. Incorrect measurement for the separation of the 
[Ne iii] has been excluded as a possible explanation, and a blending 
of the H e and the [Ne in] 3968 has been identified as the source of 
this effect. 

(viii) The measurement of A a/a using SDSS-III/BOSS spectra 
has reached the maximum precision unless better sky subtraction 
algorithms are developed. To obtain better constraints (< 10~ 6 ) us¬ 
ing the emission-line method, high-resolution spectroscopy (R ~ 
100000) is mandatory. 

(ix) We note that future large galaxies survey like eBOSS or 
DESI could provide quite stringent constraint for A a/a at low 
redshift, following our analysis of galaxy spectra taken from the 
DEEP2 survey. 
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APPENDIX A 

We publish along with this paper an electronic table with the 
combined SDSS-III/BOSS DR12 and SDSS-II/DR7 sample of 
13 175 quasars used in this work. Table A1 describes the infor¬ 
mation and format of each column. The table is available in the 
following link http://mnras.oxfordjournals.org/content/suppl/2015/ 
08/ll/stvl406.DCl/suppl _data. zip . 
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Table Al. Description of the electronic table with the combined sample (13 175 quasars) published along with the paper. 


Column 

Name 

Format 

Description 

i 

SDSS_NAME 

STRING 

SDSS-DR12 designation hhmmss.ss+ddmmss.s (J2000) 

2 

RA 

DOUBLE 

Right Ascension in decimal degrees (J2000) 

3 

DEC 

DOUBLE 

Declination in decimal degrees (J2000) 

4 

THING_ID 

INT32 

Thing_ID 

5 

PLATE 

INT32 

Spectroscopic plate number 

6 

MJD 

INT32 

Spectroscopic MJD (>55 000 SDSS-III/BOSS spectra, <55 000 SDSS-II spectra) 

7 

FIBER 

INT32 

Spectroscopic fiber number 

8 

Z_VI 

DOUBLE 

Redshift from visual inspection 

9 

Z_PIPE 

DOUBLE 

Redshift from BOSS pipeline 

10 

ERR_ZPIPE 

DOUBLE 

Error on BOSS pipeline redshift 

11 

ALPHA 

FLOAT 

A a/a from the Gaussian fits 

12 

ERR_ALPHA 

FLOAT 

Standard error for A a/a from the Gaussian fits 

13 

SN_01 

FLOAT 

S/N for the [O hi] 4960 line 

14 

SN_02 

FLOAT 

S/N for the [O hi] 5008 line 

15 

o 

1 

■n 

M 

H 

FLOAT 

Line centroid for the [O m] 4960 line 

16 

02_FIT 

FLOAT 

Line centroid for the [O m] 5008 line 

17 

ERR_01 

FLOAT 

Error on the line centroid for the [O hi] 4960 line 

18 

ERR_02 

FLOAT 

Error on the line centroid for the [O hi] 5008 line 

19 

01_AMPLITUDE 

FLOAT 

Gaussian amplitude at the centre for the [O hi] 4960 line 

20 

02_AMPLITUDE 

FLOAT 

Gaussian amplitude at the centre for the [O hi] 5008 line 

21 

01_WIDTH 

FLOAT 

Gaussian width for the [O hi] 4960 line 

22 

02_WIDTH 

FLOAT 

Gaussian width for the [O hi] 5008 line 

23 

FILE_NAME 

STRING 

File name to download from the SDSS server 





